Fusion dynamics around the Coulomb barrier 
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^ , Abstract. 

. We perform exact coupled-channels calculations, taking into account properly the effects of 

Coulomb coupling and the finite excitation energy of collective excitations in the colliding nuclei, 
^ i for three Fm formation reactions, "Cl + ^''''Bi, '^^Sc + ''^^Au, and ^'^Co + ^^^Ta. For the ^^Cl + 

■ ^"^Bi and ^^Sc + '^^Au reactions, those calculations well reproduce the experimental total fission 

cross sections, and a part of the extra-push phenomena can be explained in terms of the Coulomb 
, excitations to multi-phonon states. On the other hand, for the heaviest system, the deep-inelastic 

• collisions become much more significant, and the fission cross sections are strongly overestimated. 

j—{ I We also discuss the surprisingly large surface diffuseness parameters required to fit recent high- 

er . precision fusion data for medium-heavy systems, in connection with the fusion supression observed 

' in massive systems. 

75 ■ INTRODUCTION 

^3 : 

, The coupled-channels method has been very successful in reproducing experimental 

. ^ I cross sections for heavy-ion reactions involving medium-heavy nuclei. Particularly, in 

^ ■ many systems, it simultaneously reproduces the subbarrier enhancement of fusion cross 

sections and the shape of the fusion barrier distribution by including a few low-lying 
collective excitations of colliding nuclei and nucleon transfer channels [1]. It is now a 
standard theoretical tool to analyse experimental fusion and quasi-elastic cross sections 
at energies around the Coulomb barrier 

However, it has not yet been completely understood to what extent this method works 
for the fusion of massive systems, where the charge product of the target and projectile 
nuclei, ZpZj, is typically larger than about 1800. For those systems, other reaction 
processes, such as deep-inelastic collision and quasi-fission, come into play, and the 
reaction dynamics around the Coulomb barrier becomes much more complex than for 
medium-heavy systems [3]. In order to calculate fusion cross sections, the cornpetition 
of fusion with these other processes has to be taken into account properly Hi jsl 0], 
and the dynamics after the Coulomb barrier is overcome becomes very important. This 
is the most difficult problem in the fusion of massive nuclei, and there are still large 
ambiguities in theoretical predictions of fusion cross sections. Because of this reason, 
one often employs a simplified coupled-channels treatment for the barrier penetration 
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prior to the touching configuration |0, Rl, which is essentially based on the constant 
coupling approximation |8] or a variant 

In this contribution, we critically examine the consequences of using such a simplified 
coupled-channels framework, and point out that the exact treatment for the Coulomb 
coupling plays an important role in massive systems. We then apply the exact coupled- 
channels approach to total fission cross sections for three Fm formation reactions [10], 
where the charge product ZpZr is given by 141 1 (^^Cl -i- ^o^Bi), 1659 (^^Sc + ^^^Au), and 
1971 (^^Co -I- ^^^Ta). We will show that the small hindrance of the reduced cross sections 
for the second heaviest system compared with the lightest system can be understood 
in terms of the effect of energy loss due to the Coulomb excitation, while the fusion 
hindrance for the heaviest system exceeds that effect and an explicit treatment of deep- 
inelastic collisions is necessary. We also discuss the surface property of the nucleus- 
nucleus potential, where recent high-precision fusion data for medium-heavy systems 
systematically show that a surprisingly large diffuseness parameter is required in order 
to fit them frr. l[l^[l3i .l4]. We argue that this apparant anomaly could originate from 
the competition between fusion and the deep-inelastic processes which occur in fusion 
of heavy systems lUill . 



COUPLED-CHANNELS APPROACH TO FUSION OF MASSIVE 

SYSTEMS 

Effect of Coulomb excitations 

A hot compound nucleus formed by a fusion reaction decays either by emitting a few 
particles (mainly neutrons) and gamma rays, or by fission. The fusion cross section is 
thus a sum of the evaporation residue and the fusion-fission cross sections. Theoretically, 
it is computed as 

<Tfus(^) = ^£(2/ + i)i'fc(£,/) = ^£(2/ + i)rbp(£,/) -PcdEj), (1) 

where T\,p{E, I) is the barrier passing probability for the Coulomb barrier while Pcn{E,1) 
is the probability of compound nucleus formation after barrier penetration. 

For medium-heavy systems, the compound nucleus is formed immediately after the 
Coulomb barrier is overcome, and Pcn is essentially 1. This justifies the assumption of 
strong absorption inside the Coulomb barrier, or equivalently, the incoming wave bound- 
ary condition [2]. The coupled-channels approach has been successful here. For massive 
systems, on the other hand, PcN significantly deviates from 1, due to competition among 
several reaction processes. The fusion cross sections, therefore, appear to be hindered if 
one compares fusion cross sections with those obtained by assuming PcN = 1- Recent 
experimental data clearly indicate the strong competition between compound nucleus 
formation and quasi-fission [ 15. 16. 17]. 

In Eq. (HJ), the barrier passing probability T\,p can in principle be computed by the 
coupled-channels approach. For this, although exact coupled-channels codes are avail- 
able 10], a simplified approach has often been taken, even for fusion of massive systems 




FIGURE 1. Influence of the Coulomb couplings in fusion cross sections for the Ni + Ni (left panel) 
and the ^^Ni + ^^^Ph (right panel) systems. The double quadrupole phonon excitation in the projectile 
nucleus ^^Ni is included in the coupled-channels calculations, while the target nucleus is treated as inert. 
The solid line includes the full effect of the Coulomb excitation, while the dashed line disregards it by 
matching the numerical wave functions to the asymptotic wave functions at relatively small distance. The 
fusion cross sections without the coupling are denoted by the dotted line. 



JslSl. The simplification is achieved by using one or more of the following approxi- 
mations: i) the linear coupling approximation, where the nuclear coupling potential is 
assumed to be linear with respect to an excitation operator for intrinsic motions, ii) 
multi-phonon excitations are neglected, and iii) the eigenchannel approximation, where 
the intrinsic excitation energies are treated approximately. Since the effective coupling 
strength is approximately proportional to ZpZr for a given value of deformation param- 
eter [18], the shortcomings of these approximations become more severe for heavier 
systems. 

The validity of the first and the second approximations has been examined in detail 
in Refs. M M 

\20A . We therefore discuss the third point here. The eigenchannel 
approach is intimately related to the concept of barrier distribution [1, 8, 21]. Although 
this approach is exact when the intrinsic excitation energy vanishes, it also works well 
even with a finite excitation energy as long as the coupling potential is localized inside 
the uncoupled barrier f23\. In realistic cases, however, the coupling potential extends 
outside the barrier due to the long range Coulomb interaction. If, prior to reaching the 
barrier, there is appreciable Coulomb excitation to a collective state whose excitation 
energy is not small, this results in a decrease of the relative energy, leading to the 
reduction of fusion cross sections and a modification of effective eigen-barriers. Since 
the eigenchannel approach treats the excitation energy approximately, this effect will be 
missed in a calculation. 
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FIGURE 2. Reduced total fission cross sections for four reactions forming Fm, ^"^Cl + 2"9Bi (the 
filled circles and triangles), '^^Sc + ^^^Au (the open circles), and ^^Co + '^'Ta (the open diamonds), as a 
function of energy measured from the Coulomb barrier height. 



In order to demonstrate the effect of Coulomb excitations prior to the barrier, figure 
1 compares fusion cross sections for the ^^Ni + ^^Ni system with those for the ^^Ni + 
208p|5 system. We assume that the target nucleus is inert, and include only the double 
quadrupole phonon excitation in the projectile ^^Ni nucleus. The excitation energy 
for the single phonon state is 1.45 MeV, and we assume a simple harmonic oscillator 
coupling for the double phonon excitation. The solid and dashed lines are obtained by 
integrating the coupled-channels equations from inside the Coulomb barrier to ^max= 

50 

fm and 15 fm, respectively. The latter calculation, therefore, effectively disregards the 
effect of Coulomb excitation outside the Coulomb barrier. For a comparison, we also 
show the fusion cross sections in the no coupling limit by the dotted line. One clearly 
sees that the Coulomb coupling considerably alters the fusion cross sections when the 
target is heavy, while the difference is small for the lighter system. One also notices 
that the fusion cross sections are even smaller than the no coupling calculations for the 
heavier system at energies above the Coulomb barrier. Evidently, the Coulomb excitation 
provides another mechanism of fusion inhibition in massive systems. 



Total fission cross sections for Fm formation reactions 

Let us now discuss fission cross sections for Z=100 (Fm) formation reactions mea- 
sured by Ohtsuki et al. ifiol] . The experimental cross sections were obtained by selecting 
the fission events in the TOF vs AE plot, and thus contain both the fusion-fission and 
the quasi-fission (if any) cross sections. For the heaviest system, the fission events were 
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FIGURE 3. Comparison of the coupled-channels calculations with the experimental data for total 
fission cross sections for the Fm formation reactions. The solid line is the result of the coupled-channels 
calculations, while the dashed line is obtained without any coupling. The experimental data are denoted 
by the filled circles. 



not well separated from the deep-inelastic collision (DIG) events. For this system, the 
fission cross sections were estimated by choosing the same mass region of fission frag- 
ments as that observed in the 35,37(-i| _,_ 209gj reaction. Therefore, the experimental fission 
cross sections might be underestimated for the heaviest system if the fission events ex- 
tend towards the DIG region. Fig. 2 shows the observed reduced fission cross sections 
as a function of the difference between the center of mass energy and the average bar- 
rier energy for four systems, 35'37q + 2095^ (Z/>Zr=1411), ^^Sc + ^'^^Au (2^2^=1659), 
and ^^Go -i- ^^^Ta (ZpZ7'=1971). Notice that the reduced cross sections for the second 
heaviest system (^^Sc -i- ^^^Au) are somewhat smaller than those for the lightest sys- 
tem (^^'-^^Gl -I- ^*^^Bi), and the cross sections for the heaviest system (^^Go -i- ^^^Ta) are 
substantially hindered compared with the other two systems. 

In order to investigate whether these features can be understood in terms of the 
Goulomb excitation discussed in the previous subsection, we perform exact coupled- 
channels calculations. For this purpose, we use an extended version of the computer 
code GGFULL where the effect of finite ground state spin of colliding nuclei is 
incorporated within the isocentrifugal approximation. For the ^^Gl -i- ^*^^Bi system, we 
include three vibrational states, 5/2+ (3.09 MeV), 7/2" (3.1 MeV), and 9/2" (4.01 MeV) 



in -^^ CI as well as double octupole vibrations in ^Bi. The deformation parameters are 
estimated from the experimental B(E2) and B(E3) values. We introduce a single effective 
channel for the seven octupole states which have a character of hg ^*^^Pb(3^) in ^*^^Bi, 

and consider a harmonic oscillator coupling for the double phonon state. For the ^^Sc + 
^^^Au system, we include 5 quadrupole states which have a character of /7/2® ^^Ca(2+) 

in 45 Sc [3/2- (0.38 MeV), 5/2" (0.72 MeV), 11/2" (1.24 MeV), 7/2" (1.41 MeV), and 
9/2^ (1 .66 MeV)], and treat the target nucleus ^^^Au as a classical rotor with /32 = —0. 13 
and J84 = —0.03. For ^^Sc, since the excitation energies for the quadrupole states are 
not close to each other, we do not introduce an effective single channel, but treat them 
exactly. For the ^^Co + ^^^Ta system, we include 5 quadrupole states which have a 
character of (/7/2)"^® ^°Ni(2+) in ^'^Co [3/2" (1.1 MeV), 9/2" (1.19 MeV), 11/2" 
(1.46 MeV), 5/2- (1.48 MeV), and 7/2" (1.74 MeV)], and treat the target nucleus ^^^Ta 
as a classical rotor with 182 = 0.354 and ji^ = —0.05. More details of the calculations 
will be given elsewhere [23]. 

The results of those calculations are shown in Fig. 3. One finds that the coupled- 
channels calculations well reproduce the experimental data for the two lightest systems, 
■^^Cl + ^''^Bi and ^^Sc + ^^^Au. Especially, the reduction of cross sections in the lat- 
ter system is reproduced nicely. As we discussed in the previous subsection, this small 
hindrance of cross sections is caused by the strong Coulomb coupling to the collective 
states outside the Coulomb barrier. In contrast, the coupled-channels calculation con- 
siderably overestimates fission cross sections for the heaviest system, ^^Co -1- ^^^Ta. We 
will discuss this point in the next subsection. 



Role of deep -inelastic collision 

In his review article, Reisdorf argued [3] that cross sections for the sum of fusion and 
other damped reactions may be interpreted as the total barrier passing cross sections. 
Indeed, he found that the sum of fusion and deep-inelastic collision (DIC) cross sections 
for the ^^Ni -1- ^^"^Sn reaction could be well reproduced by the standard potential model 
(see Fig. 4). More recently, Esbensen et al. followed a similar idea and reproduced the 
experimental cross sections for the sum of fusion and DIC reactions for the same system 
with the coupled-channels approach 1I25I1 . In the semiclassical picture, deep inelastic 
collisions correspond to those trajectories which overcome the barrier in the entrance 
channel but eventually escape after appreciable interaction with the target. All of those 
considerations immediately lead to the idea that the total barrier passing cross sections, 
which the coupled-channels approach yields, may have to be compared with a sum of 
fusion, quasi-fission, and DIC, i.e., 

abp (E) = cTfus (E) + aqf (£) + ctdic (E) . (2) 

The large reduction of total fission cross sections for the ^^Co -1- ^^^Ta system, therefore, 
may be attributed to the competition between fusion, quasi-fission and DIC reactions. 

The present framework of the coupled-channels method could be used to obtain inclu- 
sive cross sections of fusion and DIC reactions, but it would be difficult to obtain them 
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FIGURE 4. Fusion cross sections (the open circles) and the sum of fusion and deep inelastic collision 
cross sections (the filled circles) for the ^^Ni + '^^Sn reaction. The solid line is a prediction of the barrier 
pen etration model with a global nucleus-nucleus potential. The experimental data are taken from Ref. 



separately. One possible way to isolate fusion from DIC theoretically is to introduce an- 
gular momentum truncation in Eq. ([T]) and exclude higher partial wave contributions, as 
was done by Zagrebaev et al. [5]. However, a large ambiguity exists in this prescription, 
since there is no clear and unique way to introduce the angular momentum truncation, 
especially at energies below the barrier ^2^]. In recent publications, Abe et al. combined 
the surface friction model with the Langevin approach in order to take into account the 
competition between fusion and DIC in the approaching phase for a formation reaction 
of superheavy elements This approach may be promising, but is essentially classi- 
cal. Computation of fusion cross sections with a quantum mechanical model under the 
influence of DIC process is still an open problem. 



SURFACE DIFFUSENESS ANOMALY IN FUSION POTENTIALS 

Let us now discuss the second subject, that is the surface property of the nucleus-nucleus 
potential and its anomaly, recently recognised in fusion reactions for medium-heavy 
systems. For scattering processes, it seems well accepted that the surface diffuseness 
parameter a should be around 0.63 fm if the nuclear potential is parametrised by a 
Woods-Saxon (WS) form. In marked contrast, recent high-precision fusion data sugges t 
that a much larger diffuseness, between 0.8 and 1 .4 fm, is needed to fit the data pH[l2il . 
This is not just for particular systems but seems to be a rather general result [14]. 

We illustrate this problem in Fig. 5 by comparing experimental data for the ^^O -i- 
208p|^ fusion reaction with various coupled-channels calculations with the WS potential. 
We include the double octupole phonon and the single 5 phonon excitations in ^*^^Pb. 
At energies well above the barrier, where the fusion cross sections Cfus is relatively in- 




FIGURE 5. Comparison of coupled-channels calculations with the experimental data for the O + 
208 fusion reaction. The double octupole phonon as well as the single 5^ phonon excitations are 
included. The solid and the dashed lines are obtained by setting the surface diffuseness parameter of 
Woods-Saxon potential to be fl=1.16 and a= 0.65 fm, respectively. The experimental data are taken from 
Ref. L2SJ- 



sensitive to the couplings, a WS potential with a diffuseness a =0.65 fm significantly 
overestimates fusion cross sections (dashed line). Changing the depth and radius param- 
eters in the WS potential is not helpful, since it merely leads to an energy shift in the 
calculated fusion cross sections without significantly changing the energy dependence. 
On the other hand, a potential with a=1.16 fm (ful l line) fits the data well. A similar 
problem also exists at energies below the barrier 11121 [l3ll. but we do not discuss it in this 
contribution. 

Up to now, several possible reasons for this anomaly have been considered. These 
include the departure of the nuclear potential from the WS form [iT', 12"], dissipation 
effects [12], and unrecognised systematic errors in experimental data [14], but none of 
them is conclusive yet. We would like to propose here another possible effect, that is, 
the competition between fusion and deep inelastic collision [|14l1 . This is motivated by 
an apparent similarity between fusion inhibition in massive systems discussed in the 
previous section (see fig. 4) and the overestimate of fusion cross sections shown in fig. 
5. To this end, let us introduce a suppression factor S which is defined as a ratio between 
the experimental fusion cross sections and the prediction of the potential model with the 
standard value for the diffuseness parameter, a ^0.63 fm. Provided that S is independent 
of bombarding energy, it can be stated that the fusion cross sections are hindered by a 
factor S for whatever reason. 

The values of S obtained for a large number of systems are shown in fig. 6 [TT]. Also 
shown are two points (large filled stars) at Z1Z2 = 1400 derived from the data for the 
^^Ni -I- ^^^'^^^Sn reactions at energies around the fusion barrier. For these systems, as we 
indicate in Fig. 4, a very significant contribution from deep-inelastic scattering has been 
observed experimentally even at energies below the barrier [24]. In this case, the value 
of S has been taken as afus divided by the sum of CTfus and ctdic- These two points do not 
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FIGURE 6. Calculated suppression factors S for fusion cross sections with respect to the prediction 
of potential model with a «0.63 fm. The large filled stars at ZiZ2=1400 refer to the reactions ^^Ni + 
Ii2,i24|gjj^ where the deep-inelastic collisions have been observed experimentally. 



seem inconsistent with the other points with large Z1Z2 in Fig. 6, which are derived from 
the fusion data only. Strongly damped reactions similar to DIC are reported to occur in 
lighter systems, for example, ^^S + ^'^Ni li29ll (ZiZ2=448) at energies well above the 
barrier. It would be interesting to see if they also occur at energies closer to the barrier. 
If they do, the values of S would be reduced below unity and might result in at least a 
partial explanation of the experimental values for a being much larger than the standard 
value, 0.63 fm. 



SUMMARY 

Extensive efforts have been made both experimentally and theoretically to understand 
fusion of massive systems, especially for synthesis of superheavy elements (SHE), 
but the reaction mechanism has not yet been completely understood. Here, we have 
performed coupled-channels calculations and pointed out that the exact treatment of 
Coulomb excitation becomes important for massive systems. We have applied the 
coupled-channels framework to three Fm formation reactions, and have shown that 
Coulomb excitation indeed provides an important mechanism of inhibition of fusion 
cross sections for transitional systems ("pre-SHE" systems) between medium-heavy 
and SHE regions. For SHE systems, we have argued that the deep-inelastic collision 
needs to be taken into account explicitly in theoretical models. However, it is still an 
open problem of how to incorporate the competition between fusion and DIC quantum 



mechanically, and further developments will be required. We have also discussed the 
large surface diffuseness problem observed in the recent high precision measurements 
of fusion cross sections for medium-heavy systems in connection with the fusion hin- 
drance in massive systems. We have argued that the competition between fusion and 
deep-inelastic collision may provide a promising avenue to explain this anomaly. In this 
connection, it would be interesting to study, both experimentally and theoretically, light 
systems near to the fusion barrier to see whether fusion was inhibited by the presence of 
Die. Especially, theoretical calculations involving friction for light systems would be of 
great interest, since experimental measurements are likely to be difficult when the fusion 
suppression factor appoaches unity. 
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